SCIENTIFIC 

REPORTS 




SUBJECT AREAS: 
BIOMATERIALS 
HIGH-THROUGHPUT SCREENING 
BIOMIMETICS 
TISSUE ENGINEERING 



open A combinatorial cell-laden gel 

microarray for inducing osteogenic 
differentiation of human mesenchyma 
stem cells 



Received 

28 October 2013 

Accepted 
8 January 2014 

Published 

29 January 2014 



Correspondence and 
requests for materials 
should be addressed to 
A.K. (alik@rics.bwh. 
harvard.edu) 

* Current address: 
Department of 
Biomedical 
Engineering, Texas 
A&M University, 
College Station, Texas 
77843. 



Alireza Dolatshahi-Pirouz 1 ' 2 ' 3 ' 4 , Mehdi Nikkhah 23 , Akhilesh K. Gaharwar 1 ' 2 ' 3 ' 5 *, Basma Hashmi' A7 , 
Enrico Guermani 2 3 , Hamed Aliabadi 2,3 , Gulden Camci-Unal 2,3 , Thomas Ferrante 1 , Morten Foss 4 , 
Donald E. Ingber w & Ali Khademhosseini 1 ' 2 ' 3 

1 Wyss Institute for Biologically Inspired Engineering, Harvard University, Boston, MA 02 1 1 5, USA, 2 Center for Biomedical 
Engineering, Department of Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA 021 39, USA, 
3 Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, MA 021 39, USA, 
interdisciplinary Nanoscience Center (iNANO), Aarhus University, Aarhus, Denmark, 5 David H. Koch Institute of Integrative 
Cancer Research, Massachusetts Institute of Technology, Cambridge, MA 02 1 39, USA, 6 School of Engineering and Applied 
Sciences, Harvard University, Cambridge, MA 021 38, USA, 7 Vascular Biology Program, Departments of Pathology and Surgery, 
Children's Hospital Boston, Harvard Medical School, Boston, MA 021 1 5, USA. 

Development of three dimensional (3D) microenvironments that direct stem cell differentiation into 
functional cell types remains a major challenge in the field of regenerative medicine. Here, we describe a new 
platform to address this challenge by utilizing a robotic microarray spotter for testing stem cell fates inside 
various miniaturized cell-laden gels in a systematic manner. To demonstrate the feasibility of our platform, 
we evaluated the osteogenic differentiation of human mesenchymal stem cells (hMSCs) within 
combinatorial 3D niches. We were able to identify specific combinations, that enhanced the expression of 
osteogenic markers. Notably, these 'hit' combinations directed hMSCs to form mineralized tissue when 
conditions were translated to 3D macroscale hydrogels, indicating that the miniaturization of the 
experimental system did not alter stem cell fate. Overall, our findings confirmed that the 3D cell-laden gel 
microarray can be used for screening of different conditions in a rapid, cost-effective, and multiplexed 
manner for a broad range of tissue engineering applications. 

The differentiation of stem cells into specialized cell types is governed by microenvironmental cues from the 
surrounding extracellular matrix (ECM) 1 " 3 , soluble factors 1 , matrix stiffness 23 , substrate topography 4,5 and 
direct cell-cell contact. These components commonly act in a synergistic manner to regulate stem cell fate 
and promote the formation of functional tissues. In particular, the three-dimensional (3D) nature of the ECM 
plays a crucial role in regulating cell behavior 6,7 . In fact, many studies have confirmed that cellular functions 
substantially deviate on 2D substrates compared to 3D microenvironments 811 . Therefore, it is important to 
develop methods of studying stem cell responses in 3D settings while controlling the presentation of other 
microenvironmental signals. An important step to meet these requirements, is the development of 3D combin- 
atorial platforms to simultaneously study stem cell differentiation in response to various cues 12 . 

Multiwell-based assays have been well accepted for screening of stem cell fate inside combinatorial environ- 
ments 7,1315 . These, platforms are typically generated by adding stem cell encapsulated hydrogels and ECM 
proteins into standard multiwell plates 13 . Other approaches have relied on the deposition and freeze drying of 
polymeric solutions into multiwell plates to generate scaffold libraries capable of screening cell-matrix interac- 
tions within 3D microenvironments 14,15 . Despite, the wide application of multiwell-based combinatorial plat- 
forms, such assays face throughput limitations due to the scarce supply of primary cells and biological signals as 
well as high reagent costs 12 . 

Recent advancements in robotic microarray technologies have enabled the development of versatile and cost- 
efficient platforms that can address the limitations of conventional screening assays 16 20 . These technologies have 
significantly contributed to our understanding of cell adhesion, proliferation and differentiation on 2D sur- 
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faces 21 " 29 . Although the 2D microarray systems have provided valu- 
able insight regarding the synergetic effects of ECM proteins on stem 
cell differentiation 21,22 ' 28 , they do not precisely mimic the in vivo 
tissue architecture. Microarray technology is currently being utilized 
to fabricate 3D miniaturized cellular platforms for drug discovery 
and toxicology research 3031 with limited focus directed toward stem 
cell differentiation 32 . For instance, a 3D cellular microgel array was 
previously created to study the effects of fibroblast growth factor-4 
and tretinoin on embryonic stem cell pluripotency 32 . Due to the few 
tested conditions, this approach did not embrace the multiplexed 
screening potential of the microarray technology. Therefore, the 
development of miniaturized platforms that enables the analysis of 
stem cell differentiation within 3D combinatorial microenviron- 
ments still needs to be fully explored 12 . 

In this work, we present a 3D cell-laden gel microarray platform 
for combinatorial screening of human mesenchymal stem cells 
(hMSCs) differentiation in response to multiple ECM and growth 
factors components. An automated printing strategy, which utilizes 
1000-fold less materials and cells, compared to conventional multi- 
well-based assays was employed to generate arrays of miniaturized 
cell-laden hydrogel constructs. Each microgel unit, composed of 
methacrylated gelatin (GE), contained living hMSCs along with 
ECM proteins and was exposed to osteogenic bone morphogenic 
proteins (BMPs). From the microarray analysis, we identified ECM 
combinations, which induced a 2-fold increase in Alkaline 
Phosphatase (ALP) expression. Furthermore, we evaluated the rel- 
evance of our platform within macroscale settings to investigate its 
translational potential. By utilizing our 3D microarray platforms, it is 
possible to efficiently screen ECM and growth factor combinations, 
which promote stem cell differentiation. We envision that, our cell- 
laden gel microarray platform could potentially accelerate the 
development of innovative and biomimetic materials for a wide 
range of applications from tissue engineering to stem cell bioengi- 
neering and drug screening. 

Results 

Fabrication of the cell-laden hydrogel microarray platform. The 

3D microarray platform was fabricated with a robotic microarray 
spotter equipped with four pins, enabling rapid printing of the gel 
spots (400 spots) in a few minutes. The schematic for the microarray 
fabrication process is shown in Figure la. Primarily, prepolymer 
solutions consisting of GE hydrogels, hMSCs and different 
combinations of ECM proteins were premixed in a 384 well plate. 
Then, the robotic microspotting system was used to print the 
prepolymer solutions on 3-(trimethoxysilyl) propyl mefhacrylate 
(TMSPMA) treated glass slides (25 mm X 75 mm) providing 
long-term surface stability of the deposited hydrogel droplets 
(Figure la). Subsequently, the deposited droplets were exposed to 
ultraviolet (UV) light to develop the miniaturized cell-laden 
constructs (Suppl. Figure la, b). Fluorescence imaging revealed a 
uniform distribution of hMSCs inside 75 um thick hydrogel 
constructs (Suppl. Figure lb), with high cellular viability after day 
1 (i.e. 91 ± 6.2% and 78 ± 6.3%) and a slightly lower viability after 7 
days of culture (Suppl. Figure 2a, b). The lower viability at day 7 is 
mainly attributed to the high crosslinking density of the 3D hydrogel 
matrix. 

To investigate the effect of ECM components on hMSCs behavior, 
a range of proteins, including fibronectin (FN, 40 ug/ml), laminin 
(LN, 40 |J.g/ml) and osteocalcin (OCN, 20 ug/ml & 40 ug/ml) were 
entrapped within the GE hydrogels (Figure b). The selected proteins 
are well known for their roles in osteogenic differentiation of 
hMSCs 33 " 38 and their affinity towards collagen I and IV 39 41 , from 
which GE is derived. We also examined the effects of growth factors, 
using bone morphogenic proteins, (BMP2 and & BMP5, 50 ng/ml), 
which promote osteogenesis via the hedgehog pathway 42 . Fluores- 
cent labeling (FN-FITC, LN-Rhodamine, and OCN-Cy5) confirmed 



the presence of FN, LN and OCN in different combinations of the GE 
microgels tested after 24 hours of encapsulation (Figure lc). 
Subsequent viability analysis demonstrated high hMSCs survival 
(i.e. 80-88%) within miniaturized 3D hydrogel constructs exposed 
to BMP2 and BMP5 (Figure Id). Overall, our preliminary results 
suggested that hMSCs can be uniformly cultured with a high cell 
survival inside the spotted 3D combinatorial microenvironments. 

Osteogenic potential of hMSCs within 3D microarrayed hydrogels. 

The osteogenic differentiation of hMSCs was evaluated by analyzing 
the expression of alkaline phosphatase (ALP) in normal (control) and 
osteogenic media using an automated high content imaging technique 
(Figure 2a). The images shown in Figure 2a suggested that combi- 
nations consisting of multiple ECM proteins lead to the highest ALP 
expression (Figure 2a). In particular, it was evident that the GE-LN- 
FN-OCN gels resulted in the most pronounced ALP up-regulation in 
both control and inducing media compared to the GE condition. 
Furthermore, the addition of BMP2 and BMP5 had less impact on 
the ALP expression levels compared to combinatorial ECM proteins 
(Figure 2a, b). 

Quantitative ALP expression analysis, depicted in a color coded 
diagram, showed that single protein conditions did not significantly 
promote the differentiation of hMSCs (Figure 2c). On the other 
hand, hydrogel constructs consisting of combinatorial ECM proteins 
(GE-FN-OCN and GE-LN-OCN-FN) induced a 2 to 4-fold higher 
ALP coverage than control condition (GE) depending on the specific 
media formulation (p < 0.05). Moreover, increasing the concentra- 
tion of OCN (20 ug/ml to 40 ug/ml) did not have a significant influ- 
ence on ALP levels. 

A one-way ANOVA analysis was performed to statistically com- 
pare all available groups for major effects on osteogenic differenti- 
ation of hMSCs (Figure 2d). This analysis indicated that out of 96 
microenvironments, 20 resulted in a major positive effect on ALP 
expression (p < 0.01). The majority of these microenvironments 
contained more than one ECM protein, thus confirming the ALP 
expression trends observed from the optical images (Figure 2a, b) 
and the color diagram (Figure 2c). Furthermore, the statistical ana- 
lysis indicated that nearly 80% of the conditions with major positive 
effect on ALP expression, contained BMP2 and BMP5 (p < 0.01). 
This further indicates that although BMP2 and BMP5 can alter the 
osteogenic fate of hMSCs, their effect on up-regulation of ALP is less 
than combinatorial ECM microenvironments. 

Osteogenic potential of specific combinations within macroscale 
hydrogels. A key aspect of our 3D cell-laden gel microarray platform 
is to predict the cell responses in larger constructs (centimeter scale) 
that can be used for regenerative medicine applications. To this end, 
specific gel compositions were selected to validate the results of the 
microarray experiments in macroscale gels. In particular, we chose to 
validate the following microarray results: a) No significant ALP 
difference was observed between GE-LN and GE gels, b) GE-LN- 
FN-OCN gels led to the highest ALP expression, and C) the addition 
of BMP2 to GE-LN-FN-OCN gels did not have a major impact on 
ALP expression. This was accomplished by examining the osteogenic 
potential of the combinations GE, GE-LN, GE-LN-OCN-FN, 
GE-LN-OCN-FN + BMP2 in control and inducing media using 
conventional macroscale encapsulation approaches. 

Optical imaging indicated an up-regulation of ALP expression in 
GE-FN-LN-OCN constructs compared to GE ones in both media 
formulations. Alternatively, in control media, the addition of LN to 
GE and BMP2 to GE-FN-LN-OCN hydrogels did not change the 
ALP expression (Figure 3a). These trends were further validated 
from high-magnification images of the hydrogel constructs as 
demonstrated in supplementary Figure 3. 

Osteogenic differentiation of hMSCs can also be identified by 
production of bone-related extracellular proteins such as osteopon- 
tin (OPN) 43 . After 14 days of culture, the GE-LN-OCN-FN and 
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Figure 1 | Fabrication of 3D hMSC-laden gel microarray. (a) A robotic microarray spotter was used to rapidly print droplets consisting of hMSCs, 
gelatin methacrylate (GE)-based prepolymer solution and various ECM proteins on TMSPMA functionalized glass slide. The printing step was followed 
by a 15 sec UV light exposure to form the miniaturized cell-laden constructs. Following printing, cell-laden gel microarrays were placed inside sealed 
chambers (Illustration made by Jeffrey Aarons). (b) Various combinations of ECM proteins and media formulations were used to conduct the 
microarrays experiments. The concentration of LN and FN was selected to be 40 ug/ml while OCN was printed at two concentrations of 20 ug/ml and 
40 ug/ml. (c) Fluorescence images of the encapsulated proteins within the hydrogel constructs after 24 hours in solution, (d) hMSCs viability within 48 
combinatorial 3D microenvironments in normal (control) media after 7 days of culture along with color-diagram displaying the quantified cell viability 
(n = 3-9). 



GE-LN-OCN-FN + BMP2 constructs showed a higher OPN cov- 
erage compared to the GE gels in both media formulations 
(Figure 3b). Consistent with these observations, the quantification 
of ALP activity and OPN coverage demonstrated that the GE-LN- 
FN-OCN gels led to a significantly higher ALP expression, while 
GE-LN-OCN-FN + BMP2 in control media up-regulated OPN 
expression (Figure 3d). Notably, the ALP activity observed in macro- 
scale setting (Figure 3c) was closely correlated with the microarray 
results as confirmed by a linear regression fit (Suppl. Figure 4). 

Another hallmark of osteogenic differentiation is the production 
of calcified matrix 44 . Therefore, we conducted further analysis to 
identify the ECM combinations, which promoted matrix calcifica- 
tion. As expected, there was no notable calcification in the tested 



conditions under control medium when compared to the inducing 
medium, as detected by Alizarin Red staining (Figure 4a). However, 
the production of calcified matrix in GE-LN-FN-OCN was signifi- 
cantly higher compared to GE hydrogels in inducing media after 14 
and 25 days of culture. After 25 days of culture, the GE-LN-FN-OCN 
constructs displayed dense microstructures that were heavily stained 
by Alizarin Red in accordance with formation of bone -like calcified 
structures (Figure 4a). Theses observations were further confirmed 
through quantification of Alizarin Red stain as shown in Figure 4b. 

Earlier studies have demonstrated that production of mineralized 
matrix also induce a transformation from transparent into trans- 
lucent/opaque hydrogel structures 45 . Consistently, the hydrogels cul- 
tured in non-inducing media were nearly transparent, whereas 
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Figure 2 | ALP expression within the 3D hMSC-laden gel microarray platform, (a) Optical images of ALP expression of hMSCs within various 3D 
combinatorial microenvironments. (b) High magnification images indicating ALP expression inside the representative constructs, (c) Color-diagram 
representing the normalized ALP coverage within each 3D microenvironment. Groups with the same color code were not statically different from one 
another (n = 9-22) (d) The scaled estimates of major effects on ALP expression levels calculated using one-way ANOVA analysis. 



hydrogels in inducing media displayed opaque-white microscale fea- 
tures after 25 days in culture (Suppl. Figure 5). In particular, GE-LN- 
FN-OCN gels exhibited highest level of white coloration (Suppl. 
Figure 5) caused by clusters of micrometer- sized nodules, similar 
to the calcified microparticles seen in Figure 4a. 



To investigate the mineral phase of the deposited matrix, a series of 
Raman spectra was acquired from the macroscale constructs. As 
shown in Figure 4c, a sharper and more intense hydroxyapatite 
(HA) peak was detected in GE-LN-FN-OCN compared to GE gels 
in inducing media confirming the formation of crystalline HA 46 . 
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Figure 3 | Microarray hit combinations stimulates early osteogenic differentiation of hMSCs within the macroscale GE gels, (a) ALP expression after 7 
and 14 days of culture, (b) Osteopontin (OPN) expression and DAPI staining of cell nuclei along with (c) quantified ALP activity after 7 and 14 
days of culture (n = 7-15). (d) Quantification of the normalized OPN coverage after 14 days of culture (n = 9-15). 



Alternatively, there was no HA peak present in the spectral analysis 
conducted on the control conditions. Matrix mineralization is well 
known to enhance the mechanical properties of hydrogels 45 . 
Therefore, an unconfined compression test was performed to evalu- 
ate the mechanical properties of our engineered tissue constructs 
after 25 days of culture. The results indicated a significant increase 
in the compressive modulus of GE-LN-FN-OCN gels (16.3 kPA) 
compared to the GE gels (4.8 kPA) in inducing mediuma 
(Figure 4d). In addition, all the constructs tested in the inducing 
media exhibited a significantly higher compressive modulus com- 
pared to the samples cultured in control media. In summary, our 
study confirmed that GE-LN-FN-OCN stimulates osteogenic differ- 
entiation of hMCSs and promotes bone mineralization as predicted 
by the 3D cell-laden gel microarray platform. 



Discussion 

Fabrication of a miniaturized cell-laden gel microarray platform 
requires an optimal printing technique, appropriate biomaterials 
(i.e. hydrogels) and functionalized substrates to ensure high cell sur- 
vival and experimental stability for long-term studies. We utilized a 
robotic microspotter system equipped with multiple pins for rapid 
deposition of the cell-laden hydrogel constructs inside a temperature 
and humidity controlled chamber to keep the cells viable during the 
fabrication process. Furthermore, TMSPMA functionalized substrates 
and GE prepolymer solution were used to ensure high cellular viab- 
ility, successful ECM protein and cellular encapsulation and platform 
stability up to 7 days of culture. We noted that it was challenging to 
maintain the stability of our microarray platform beyond 7 days of 
culture due to hydrogel degradation most likely caused by hMSC 
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Figure 4 | Microarray hit combinations stimulates macroscale bone-mineralization, (a) Formation of calcified matrix by hMSCs investigated using 
Alizarin Red S staining after 14 and 25 days of culture; white circles denotes regions of intensified Alizarin Red S staining, (b) Quantified amount of 
calcified matrix evaluated by Alizarin Red S content inside the hydrogels using a colorimetric assay (n = 4-8). (c) Raman spectra of the hydrogel 
constructs after 25 days of culture, (d) Unconfined mechanical testing to evaluate the compressive modulus (kPa) of the macroscale hydrogels upon 
formation of mineralized matrix (n = 4). 



spreading, migration and collagenase secretion of hMSCs. We could 
potentially improve this by adding naturally derived hydrogels with 
low degradation rate (i.e. alginate, hyaluronic acid) to the prepolymer 
solution. 

To date, significant efforts have been devoted to develop 2D 
microarray platforms. However, by using 2D settings, it is difficult 
to fully understand the role of matrix cues on stem cell differentiation 
inside native-like 3D microenvironments 6 . Moreover multiwell- 
based arrays do not enable high-throughput studies. Therefore, we 
believe that the 3D microarray platform presented here better simu- 
lates in vivo responses of stem cells, while allowing cost-efficient and 



high-throughput analysis. In fact, there is growing evidence that stem 
cells elucidate more biomimetic behavior inside 3D microenviron- 
ments compared to 2D surfaces 811 . Specific examples include osteo- 
genic 8 , and hepatic differentiation 9 of embryonic stem cells (ESCs) 
within 3D microenvironments. 

Our interest in bone tissue engineering prompted us to focus on 
ECM proteins and growth factors that are well known to induce 
osteogenic differentiation of hMSCs. Both FN and LN promote 
osteogenic differentiation through binding to integrins a5FJl and 
oc3pT located on the membrane of hMSCs 34,35 . Similarly, OCN plays 
an important role in the formation of mineralized matrix through 
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modulation of HA crystal morphology and growth 33 . Our findings 
indicated that, combinations of FN, OCN and LN resulted in a higher 
ALP expression compared to individual incorporation of ECM pro- 
teins (Figure 2). It is widely accepted that functional properties of 
ECM proteins can be changed through structural alterations directed 
by protein-protein interactions 39,47,48 . Such structural alterations 
could ultimately expose hidden osteogenic regions on FN and LN, 
thereby enhancing the osteogenic differentiation of hMSCs. 
Therefore, it is reasonable to assume that these synergistic interac- 
tions may be responsible for the up-regulation of ALP observed in the 
microarray experiments. Similar synergistic interactions were also 
observed by Flaim 21 and Fang 13 et al, upon screening of ESCs res- 
ponses to different ECM combinations on 2D micro- and 3D multi- 
well-based arrays, respectively 21 . They reported that the crosstalk 
between collagen, LN and FN could lead to up-regulation of stem 
cell markers such as albumin (hepatic marker), osteocalcin (bone 
marker) and von Willebrand factor (endothelial marker). These 
studies reinforce the fact that ECM proteins significantly influence 
the biological functionality of stem cells through structural cues that 
are directed by protein-protein interactions. 

The applicability of our developed 3D microarray platform was 
further investigated by evaluating the osteogenic differentiation 
of hMSCs within macroscale gels containing ECM proteins. 
Significantly higher ALP expression was observed in GE-LN-FN- 
OCN compared to GE gels after 7 and 14 days of culture. The 
increase in ALP expression is known to promote bone mineralization 
by removing extracellular inorganic pyrophosphate (ePPi), a cal- 
cification inhibitor, and releasing inorganic phosphate 49 . Calcium 
quantification and Raman spectral analysis of the gel constructs 
along with mechanical tests performed on the GE-LN-FN-OCN gels 
further confirmed that 'hit' combinations retrieved from the pro- 
posed microarray platform could potentially enhance macroscale 
bone formation. In summary, our findings demonstrated that the 
proposed 3D microarray platform provide a unique setting to screen 
cell-matrix interactions in a high throughput manner for tissue 
engineering applications (i.e. bone) with potential translational 
impact. 

Conclusion 

In this work, a microarray platform consisting of miniaturized 
hMSC-laden hydrogel constructs was developed for high throughput 
and cost efficient screening of stem cell-matrix interactions. The 
microarray system allowed for rapid printing of hMSCs inside vari- 
ous 3D milieus with high cell survival rate. The utility of the micro- 
array platform was further demonstrated by evaluating the 
osteogenic potential of hMSCs inside 96 different miniaturized 3D 
niches. A range of ECM proteins were identified, which induced 
significantly higher ALP expression in combinatorial settings com- 
pared to constructs containing single protein. The hit combinations, 
identified using the microarray platform, were reinvestigated in 
macroscale setting confirming their potential in the formation of 
mineralized bone tissue. Our findings indicated that the 3D minia- 
turized stem cell niches could be used to predict cellular responses 
and explore appropriate matrix components for a wide range of 
tissue engineering applications. 

Methods 

hMSCs cell culture. Bone-marrow derived hMSCs were purchased at passage 2 from 
a commercial source (Lonza, MD) and expanded in a T-150 flask in mesenchymal 
stem cell growth media (MSCGM) (Lonza, MD) to passage 3-5 before use. Cells 
where passaged when 80-90% confluent in accordance with the manufactures 
instructions and the growth medium was changed twice a week. For osteogenic 
differentiation studies, the cells were cultured in osteogenic induction media (Lonza, 
MD) containing mesenchymal cell growth media supplemented with 
Dexamethasone, L-Glutamine, Ascorbate, Pen/Strep and |3- Glycerophosphate. 

Methacrylated gelatin synthesis. GE was synthesized according to previously 
published method 50 . Briefly, 10 g of gelatin (Sigma) was diluted in 100 mL PBS at 



50°C. Then, methacrylation was induced by slowly adding 8% (v/v) methacrylic 
anhydride to the gelatin solution and reacting at 50°C for 3 h. The final GE solution 
was dialyzed using a 12-14 kDa dialysis membrane for one week. Following the 
dialysis process, the solution was frozen at — 80°C and later lyophilized for 4 days 
before the final use for the microarray experiments. 

Microarray printing. The microarray printing process was performed on glass slides 
functionalized with 3-(trimethoxysilyl)propyl methacrylate (TMSPMA) (Sigma, 
MO). Prior to printing, different combinations of proteins and cell contained 
hydrogel solution were mixed and injected into 384- well plates. Specifically, the 
mixtures contained Phosphate Buffered Saline (PBS) (Invitrogen, NY), hMSCs, GE 
prepolymer solution and various combinations of ECM proteins including FN 
(Abeam, MA), OCN (Abeam, MA) and LN (BD Biosciences, CA). hMSCs were 
trypsinized and suspended in a solution containing 5% (w/v) GelMA and 0.5% (w/v) 
Irgacure 2959 photoinitiator ((2-hydroxy-l-(4-(hydroxyethoxy) phenyl)-2-methyl- 
1-propanone) (CIBA Chemicals, CH) to reach a final cell concentration of 8 X 
10 6 cell/ml of hydrogel solution. FN and LN concentrations were selected to be 40 p.g/ 
ml, while two different concentrations of 40 ug/ml and 20 ug/ml were used for OCN. 
The mixtures were printed on the TMSPMA coated glass slides using a SpotBot 3 
contact microarrayer with four pinheads (Arrayit, CA). During the printing process, 
the temperature and humidity inside the microarray chamber were kept above 90% 
and 22°C respectively. The high level of humidity was used to avoid polymer 
evaporation from the printed spots and potential cell death. The deposited microspots 
were subsequently crosslinked using UV light (320-500 nm, 800 mW power, 10 cm 
height) for 15 seconds. Following the printing process, the samples were cultured in 
hybridization chambers (Grace Bio, OR) containing MSCGM media. After 24 hours, 
4 different combinations of BMP2 and BMP5 (Abeam, MA) were added to each 
chamber. For osteogenic differentiation experiments, the MSCGM media was 
replaced with osteoinductive media after 24 hours of culture. We assured that the 
results obtained from the microarray platform was not influenced by cross- 
communication between the neighboring gels, by printing each microgel 
combination in isolated 3X3 spot quadrants. Moreover, the experiments conducted 
on glass slides containing scrambled microgel blocks also resulted in similar findings 
(data not shown), which may have minimized the spot-to-spot crosstalk. Images of 
the printed microarrays were acquired at 5X magnification using a Zeiss Axio 
Observer Zl microscope equipped with a programmable motorized stage for high- 
content imaging. Images analysis was performed in a high throughput fashion using 
National Institute of Health (NIH) ImageJ software (v. 1. 4). 

Confocal imaging. Confocal images were acquired using a Leica TCS SP5 microscope 
equipped with a 405 nm diode and a white light laser tuned to 495 nm, a UV 
corrected HCX PL APO CS 20.0X/0.70N.A. air objective, and Leica hybrid detectors 
collecting at 415-485 nm and 505-596 nm, respectively, with 1.51 X 1.51 X 2.9 Jim 
voxels (Leica Microsystems, IL). The 3D constructed images were acquired using 
Imaris v7.6.1 software (Bitplane Inc., CT). 

Viability analysis. A Live/Dead Kit (Invitrogen) was used to quantify the percentage 
of live cells after 1,3,5 and 7 days of culture by exposing the printed microarrays slides 
tocalceinAM(Cl) (0.5 ill/ml) and ethidium homodimer (ETD) (2 ul/ml)for30 min. 
The stained slides were imaged at 5 X magnification using an inverted fluorescence 
microscope (Zeiss Axio Observer Zl, Zeiss, USA). The obtained images were then 
analyzed using ImageJ v. 1. 4 software. 

Histology. Hydrogel samples were fixed in 4% (v/v) paraformaldehyde (PF) 
(Electron Microscopy Sciences, PA) and placed in 30% (w/v) sucrose overnight. 
Samples were then embedded in OCT and cryosectioned (Leica Cryostat) at 20 um. 
Each hydrogel section was then prepped according to the subsequent 
immunostaining procedure provided by the manufacturer. 

Osteogenic differentiation. After fixing the samples in 4% (v/v) PF for 20 min, the 
osteogenic differentiation of the cells was examined using ALP and Osteopontin 
(OPN) expression as the early differentiation markers. Alizarin Red S was further 
used to define the matrix calcification as a late marker. A standard substrate 
containing BCIP/NBT (MP biomedical, OH) was used to evaluate ALP expression. 
We gently added the solutions to each slide and incubated them for 90 min. The slides 
were then rinsed three times in PBS for subsequent imaging. The ALP activity was 
evaluated using a commercially available colorimetric Assay kit (Abeam, MA). 
Briefly, the gels were digested with 100 ug/ml dispase II (Sigma, USA) and added to a 
96-well plate and optical density (O.D.) measurement was performed at 405 nm 
using a SpectraMax M5 (Molecular devices, CA) spectrometer. From standard curves, 
we calculated the ALP activity within each representative hydrogel construct. The 
ALP activity was finally normalized to the DNA content using a Quant-iT PicoGreen 
kit (Invitrogen, NY) according to the manufacturer's protocols. For OPN staining, the 
samples were primarily blocked for 2 hours using 2.5% (v/v) goat serum (Invitrogen, 
NY), 1% (w/v) bovine serum albumin (BSA)(Sigma) and 0.3% (v/v) Triton (Sigma). 
The samples were then exposed to 1 : 200 dilution of primary polyclonal rabbit anti- 
OPN antibody (Abeam, MA) in 1% goat serum in PBS for overnight at 4 L C. 
Following, a 1 : 200 dilution of Alexa Fluor-555 conjugated goat anti-rabbit secondary 
antibody (Invitrogen, NY) in 1% goat serum in PBS was added to the slides for 2 hours 
at room temperature. Between each step the slides were gently washed three times in 
PBS. The samples were washed repeatedly for three hours in ddH 2 0 and afterwards 
incubated in ddH 2 0 overnight to remove excess stains. Quantitative analysis of 
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Alizarin Red S was performed by dissolving the stained deposits in 10% (v/v) acetic 
acid for 30 min followed by 30 min in 10% ammonium hydroxide (w/v) to ensure the 
pH-value was maintained in the range 4.1-4.3. Finally the OD 405 was measured with a 
SpectraMax M5 {Molecular devices, CA) and compared with the standard curve to 
extract the Alizarin Red S content from each sample. In all the differentiation readouts 
we used GE constructs cultured in non-inducing medium as our negative control. 

Raman spectroscopy. Mineralization of the matrix was investigated using R200-L 
SENTERRA Raman microscope equipped with Olympus BX51 microscope stand. 
The sample was exposed to He-Ne laser at 785 nm wavelength with a power of 5 mW 
as the excitation source. The final spectrum was collected through averaging of five 
samples with 8 seconds in between acquisition time. 

Mechanical testing. Mechanical properties of hydrated macroscale hydrogel 
constructs were determined by unconfmed compression testing using Instron 5943 
Materials Testing System Capacity (Norwood, MA, USA) equipped with a 50 N load 
cell. The uniaxial compression testing was performed with a crosshead speed of 
1 mm/min on circular samples with 8 mm diameter and 1.0-1.4 mm thickness until 
point of failure. The 5-15% strain region was used to measure the compressive 
modulus of the samples. 

Statistical analysis. To assess the level of significant differences between different 
experimental groups in the microarray analysis, a one-way ANOVA was performed 
followed by a post hoc test. We used a Tukey- Kramer all pairs post hoc test when the 
groups followed a normal distribution and had equal variances, otherwise a Steel- 
Dwars all pairs non- parametric post hoc test was used. Groups, which were not 
statistically different, were organized in the same color group represented in a color 
diagram. To reveal combinations, which had a major osteogenic effect, we carried out 
a scaled effect screening within various combinations using a one-way ANOVA 
analysis. 

For the macroscale experiments, we followed the same approach when comparing 
more than three groups. Otherwise either a two-tailed t-test was performed, when the 
data was normally distributed, or a non-parametric Mann- Whitney test was carried 
out. In general, averages were considered to be significantly different when the p- 
value was less than 0.05. All the data were presented as average ± standard deviation 
(SD). All the statistical analysis was carried out using JMP software (vs. 10). 
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